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An arbitrary quantum-optical process (channel) can be completely characterized by probing it 
with coherent states using the recently developed coherent-state quantum process tomography 
(QPT) [Lobino et ai, Science 322, 563 (2008)]. In general, precise QPT is possible if an infi- 
nite set of probes is available. Thus, realistic QPT of infinite-dimensional systems is approximate 
due to a finite experimentally-feasible set of coherent states and its related energy-cut-off approxima- 
tion. We show with explicit formulas that one can completely identify a quantum optical Gaussian 
process just with a few different coherent states without approximations like the energy cut-off. For 
tomography of multimode processes, our method exponentially reduces the number of different test 
states, compared with existing methods. 
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I. INTRODUCTION 

One of the basic problems of quantum physics is to 
predict the evolution of a quantum system under cer- 
tain conditions. For an isolated system with a known 
Hamiltonian, the evolution is characterized by a unitary 
operator determined by the Schrodinger equation. How- 
ever, the system may interact with its environment, and 
the total Hamiltonian of the system plus the environment 
is in general not completely known. The evolution can 
then be regarded as a "black-box process" [ll-Q which 
maps the input state into an output state. An impor- 
tant problem here is how to characterize an unknown 
process by testing the black-box with some specific input 
states, which is referred to as quantum process tomogra- 
phy (QPT) (for reviews see Refs. [1,[|). 

QPT can be understood as the tomography of a quan- 
tum channel since any physical operation describing the 
dynamics of a quantum state can be considered as a chan- 
nel Q. By contrast, the goal of quantum state tomog- 
raphy (QST) is the reconstruction of an unknown state 
(i.e., its density matrix) by a series of measurements on 
multiple copies of the state (for a review see Ref. [J]). 
Both QPT and QST are essential tools in quantum engi- 
neering and quantum information processing. 

A few methods for QPT were developed including the 
standard QPT [l|, Q], ancilla-assisted QPT JlBjp, di- 
rect characterization of quantum dynamics jlOl lll| , and 
coherent-state QPT ^^. There are dozens of propos- 
als and experimental realizations of QPT for systems 



with a few qubits. These include the estimation of quan- 
tum optical gates M, 11 , [l3l - [l9f , liquid nuclear- magnetic- 



11,1131- 

resonance gates |20l - l22i superconducting gates 
(fo r a review see Ref. [29|) and other solid-state gates (30|- 
|32| . ion-trap gates [3^, [3J], or the estimation of the 
dynamics of atoms in optical lattices [35| . By con- 
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trast, there are only a very few experimental demonstra- 
tions of QPT for infinite-dimensional systems (see, e.g., 
Ref. [11). 

Any physical process can be described by a completely 
positive map e. Such a process is fully characterized if 
the evolution of any input state pm is predictable: pout = 
e(pin). In general, QPT is very difficult to implement 
in high-dimensional spaces, and, more challengingly, in 
an infinite-dimensional space, such as a Fock space [9|, 
Il3 |. Recently, Ref. [ll showed QPT in a Fock space for 
continuous variable (CV) states. Two conclusions can be 
drawn ^^■. {i) If the output states of all coherent input 
states are known then one can predict the output state 
of any input state; (m) By taking the photon-number- 
cut-off (or energy cut-off) approximation, one can then 
characterize an unknown process with a finite number of 
different input coherent states (CSs). 

It is an interesting question to identify an exact QPT 
with finite coherent states. If the process is completely 
unknown then QPT with finite coherent states is impos- 
sible. However, if some of the constraints of the quantum 
process are known then QPT can be simplified and, thus, 
effective. In quantum optical process, Gaussian maps are 
the most common. In this article we show that if a cer- 
tain quantum process is known to be Gaussian then an 
exact QPT can be performed with only a few different 
coherent states. 

It is worth noting that there is some analogy between 



QPT and QST, especially for quantum optical Gaussian 
processes (channels) [y, jSg and Gaussian states [33, [S^ ■ 
This analogy can be seen, e.g., by comparing correlations 
between observables encoded in the covariance matrices, 
which completely describe a Gaussian object (either a 
quantum state or process). Thus, tomographies of quan- 
tum Gaussian systems arc effectively finite-dimensional 
with their covariance matrix having a physical meaning 
analogous to a finite-dimensional density matrix. 

As has been shown in Ref. [^ , QPT can be performed 
with two-mode squeezed vacuum (TMSV) for any un- 
known process. However, TMSV are not so easy to ma- 
nipulate in practice. Especially, because this involves 
quantum tomography of entangled states, which is not 
an easy task. 

Here, we show that based on existing results [j], 
by using the standard quantum-optical Husimi Q- 
rcprcsentation, one can perform QPT with only a few CSs 
without entangled ancillas for quantum optical Gaussian 
processes. The method described here has several advan- 
tages. First, it presents explicit formulas without any ap- 
proximations, such as the photon- number-cut-off approx- 
imation. Second, it requires only a few different states 
to characterize a process, rather than all CSs. Third, 
for multimode Gaussian process tomography, the num- 
ber of input CSs increases polynomially with the number 
of modes, rather than exponentially. Fourth, it uses the 
Husimi Q-functions only, which is always well-defined for 
any state without any higher order singularities in the 
calculation. 

The paper is organized as follows: In Sec. H, the above- 
mentioned QPT without ancillas is proposed for single- 
mode Gaussian processes. A simple illustrative example 
of the method is discussed in Sec. H.A. A generaliza- 
tion of our QPT for a multi-mode case is presented in 
Sec. HI. For comparison, we review the existing results 
about QPT based on entangled ancillas in Appendix A. 



II. ANCILLA-FREE GAUSSIAN PROCESS 

TOMOGRAPHY WITH A FEW COHERENT 

STATES 

Here we present the main result of this paper, which 
is an efficient tomography of Gaussian processes probed 
with only a few single-mode coherent states without the 
assistance of ancillas. 

As shown in Ref. [l^, if we only use CSs in the test, 
the tomography of an unknown process in Fock space 
requires tests with all CSs. Though this problem can 
be solved by taking the photon-number-cut-off approx- 
imation, in a quantum-optical process associated with 
intense light, one still needs a huge number of different 
CSs for the test. Here we show that the most important 
process in quantum optics, the Gaussian process |6l.[36|. 
can be exactly characterized with only a few CSs in the 
test. 

A Gaussian process maps Gaussian states into Gaus- 



sian states [39f . Therefore the Husimi Q-function of the 
operator p^ must be Gaussian: 

Qp.(z:,Z*,Z„Zb) =exp(co + L + L^ + 5 + 5t + S'o), (1) 

where 

b 
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Before testing the map, all these are unknowns. The 
normally-ordered form of the density operator p^ is 



Pe ~- Qp,(a^&^a, &) 



(2) 



corresponding to Eq. ([T|) but with variables 
{Zl,Zl,Za,Zi,) replaced by {a\b\a,b). The nor- 
mal order notation : . . . : indicates that any term inside 
it is reordered by placing the creation operator in the 
left. For example, : aba'^h^a := a'^h^a^h. 

The output state from any single-mode input coherent 
state \u) (in mode b) is 

Pu^ira[{\u*){u*\)^®I-pe]. (3) 

where p^ is given by Eq. ^. Its Husimi Q-function is 

QpAZbiZb) = Qp^{u,Zb,u*,Zb) 

= exp{cu+Lu + Ll + R + R^ +Ro), (4) 

where 

Lu = (r& + U*Xab + uYab)Zb, 

R = ZbXbbZb/'i', Ro = Z^YbbZb, 
and Cu is determined by cq, Fa, Xaa, and Yaa- Explicitly, 

C„ = Co + Re {2TaU* + U*XaaU* + uYaaU*) . (5) 

The quadratic functional terms {R,R^,Rq) in the expo- 
nent in Eq. (J4|) are independent of u; these terms must 
be the same for the output states from any input CSs. 
Therefore, these can be known by testing the map with 
one coherent state. Thus, we do not need to consider 
these terms below. 

Now suppose that we test the process with six different 
CSs, \ai), and i = 1,2,--- ,6. Assume also that the 
detected Husimi Q-function of the output states is 

Qp^, {Z;,Zb) = cxp(c, + D, + dI+R + R^ + Ro), (6) 



where Di — diZi, is the detected (hence known) linear 
term. We note that there are available efficient methods 
of Gaussian QST based on homodyne detection, which 
enable the estimation of the Wigner function or, e quiv - 
alently, the Husimi Q-function for Gaussian states |38|. 
According to Eq. (|4]), the Q-function of the output state 
from the initial state \ai) of mode b must be 



?P„. (^fc.^b) = QpA<^i^Z^,ai*,Zb). 



(7) 



Therefore, we can derive self-consistent equations by us- 
ing the detected data from p^. and setting u = a^ in 
Eq. Q: 



U - A, 



(8) 



where Li , Ca^ are just L^, c„ , respectively, after setting 
M = tti in Eqs. ©-([S]); Di and Ci are known from tests. 
Explicitly, 



Li = (Tb + ai*Xab + aiYab)Zb- 
The first part of Eq. ([5]) causes: 



K ■ (Eb, Xab, Yab) = d, 



(9) 
(10) 



where 





A' = 



There are three unknowns (E^, Xab, and Yab) with three 
equations now. Wc find 



(Efc, Xab, Yabf = K-^d 



fll) 



If the Gaussian process is known to be trace-preserving 
then Eq. (jlip completes the tomography: up to a nu- 
merical factor, we can deduce all the output states of the 
other input CSs, \ai), for i = 4, 5, 6. The term c; can be 
fixed through normalization, which is determined by the 
quadratic and linear functional terms in the exponent of 
the Q-functions. Knowing these {ci}, one can construct 
Pe completely as shown below. 

Eor any map, Ci can be known from tests with \ai), for 
i = 1, 2, • • • ,6. We then have 



•J ' 1^0, i aj t ai -^aa, -^ aa^ ^ aa) 



(12) 



where 
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for the second part of Eq. 



(co, Tq 






dHl). Thus 



C2 
C3 
C4 



(13) 



Theorem: Given K and J defined by Eqs. (|T0)) - ((T2)) then 
the QPT of any single-mode Gaussian process in Fock 
space can be performed with six input CSs, when det K ^ 
and det J 7^ 0. The QPT of any trace-preserving single- 
mode Gaussian process in Eock space can be executed 
with three input CSs, when det/\ ^ 0. 

For example, one can simply choose ai = 0, a2 = 1, 
as = i, ai = —1, as = — «, and a^ = 1 + i. One finds 

Co = ci, Tb = di, 

Ta = T (C2 + «C3 - C4 - IC^) , 

Xab = -;^[-[l + i)di+d2+id3], 
Yab = ■^[-il-i)di+d2-id3], 



(14) 



Y, 
X. 



aa - t(c2 + C3 + C4 +C5) - Ci, 



aa - - \Cl - C3 + C4 - C5 + 2i(ci - C2 - C3 + Cg)] , 



where \di\ and {ci} are defined in Eq. (jB]). 



A. Example: Output state of a beam-splitter 
process 

As a check of our conclusions, we calculate the output 
state of a beam-splitter (BS) process as shown in Fig. 1. 
The BS has input modes h and c and output modes h' 
and c'. Regarding this as a black-box process, the only 
input is mode h and the only output is mode 6'. We 
set mode c to be the vacuum. The BS transforms the 
creation operators of modes h and c by: 



C/bS (&^: 



c^)t/, 



-1 
BS 



-(^^ 



A\ 



Af, 



BS, 



(15) 



where A/bs = 



If wc test such a pro- 



cos ty sm t 
— sin 6* cos ( 

cess with a coherent state \ai), we shall find pa^ = 
jofi cos6')(Q;i cos6'|. Comparing this with Eq. ©, we have 
di = a* cos 6' and q = -|ai cos 6lp. Using Eqs. ([TT|) -([T ^ . 
we find 



Ybb = -1, ATah =cos( 
Ta = Tf, = Yab = A(jQ 



X, 



bb 



co = 0. 



(16) 



Therefore 

Pe =: exp(a^5t cos 6' — a^a cos^ 9 — b^b + abcosO) : . (17) 

With this we can predict the output state of any input 
state, for example the squeezed coherent state (squeezed 
displaced vacuum) 

\ar, Z)) = cxp lUb" - 6t2)l cxp (Z5t - Z*b) |0), (18) 



where r is real. According to Eq. (|A8|) . 

p^ = tr, [(|e(r, Z*))(e(r, Z*)\)a ® h ■ Pe] ■ (19) 



Beam-splitter 
b 




Input of the 
Gaussian Map 




Output of the 
Gaussian Map 



which is a {2k + 1) x (2fc + 1) matrix, since each a^ here 
is a fc-mode row vector. Moreover, d is here a {2k + 
1) X fc matrix as (f" — {d[ , d^, • • ■ j^ifc+i/j with di — 
{dii,di2,--- ,dik)- Similarly, 



FIG. 1: Gaussian map constructed by a beam-splitter. Here, 
we assume that the input mode c is in the vacuum. 



As a result, 

Qp^ {Z; ,Zb)^C exp{Hi -n2 + 'H3+n4), 
where C is the normalization factor, and 



(20) 



Hi 



|Z6p(tanhVsin^0-l)/5, 



•Hs = {Z^ + Z;^)tanhrcos'^0/{2g), 

Hz = ZbCose{Z* - Zt!mh.rsm^e)/{gcos\Yr), 

■Hi = Zb*cos6l(Z-Z*tanhrsin2 6')/(gcoshr), 

and (7 = 1 — tanh^ r sin^ 9. This is the same result ob- 
tained from direct calculations using Eq. ([T5|) . 



III. EFFICIENT MULTIMODE-GAUSSIAN QPT 

Multimode Gaussian QPT has many important appli- 
cations. For example, it applies to a complex linear op- 
tical circuit with BSs, squeezers, homodyne detections, 
linear losses, Gaussian noises, and so on. Consider now 
a Gaussian process acting on a /c-mode input state (in 
modes foi, 621 ■ ' ' 7 ^fc)? with outcome also a fc-mode state. 
Even though other methods [IJl can also be extended to 
the multimode case, the number of input states required 
there increases exponentially with the number of modes 
k, because the number of ket-bra operators |{"-i})({m-i}l 
in Fock space increases exponentially with fc. As shown 
below, the number of input states in our method increases 
polynomially. 

A fc-mode QPT can be tested with fc-mode CSs, if the 
process is Gaussian. The main Eqs. (|lip - p^ still hold 
after redefining the notations there. First, Fa, F;,, u, a^, 
di, Za, and Zf, arc now fc-mode vectors. For example, 

\ai) = \aii,ai2, ■ ■ ■ ,a,k), 
di = {dii,di2, ■ ■ ■ ,dik), Zb = {Zhi,Zh2, ■ ■ ■ ,Zhk), 



Xxy is now a fc X fc matrix. 



and so on. Following Eq. ([T]), 
for X ^ X OT Y with x — a,b; y ~ a,b. We still apply 
Eqs. dm-dll]) to calculate {Ff,, Xab, Yab} and {Fa, Xaa, 
Yaa}, respectively, but keep in mind that the matrices K, 
J and symbols d, c are now redefined. There are (2fc-|-l)fc 
unknowns in (Tb, Xab, Yab)- We need (2fc-|- 1) different 
CSs of fc-mode to fix these unknowns. Now we have 



/I 



K = 



1 a*n 



1 at 



ai ^ 

a2 

ctk ) 



(21) 



/I 



J = 



a.2 



\af \a\ 



\ 



2"2 



laal 



\\ a'j^ aN \ 



*2 



(22) 



2^N 2"JV 



CO. 



Wn? ) 



which is now a. N x N matrix, and N = {k + l)(2fc + 
1), since of and Ictip here are row vectors of of = 
{Eii, Ei2, • • • , Eik) and la^p = {En, Ei2, ■■■ , E^k), and 
each element of Eim (or Eim) is a vector with (fc — to+ 1) 
modes (or fc modes), as 



Ei: 



( 2 



iOik,al), 



/ * * * * \ 

[ainian , 0^ijji(Xi2, • ■ ' , Q^mi^j.fc— 1 ; aij^iCXif^j. 



Obviously, c is a column vector with N elements. There- 
fore we conclude with this: 

Corollary: Any fc-mode Gaussian QPT can be per- 
formed with (fc + l)(2fc + 1) different CSs of fc-mode; 
or with (2fc -|- 1) different CSs of fc-mode if the process is 
trace preserving. 



IV. CONCLUSIONS 

In summary, we have presented explicit formulas for 
the tomography of quantum optical Gaussian processes 
probed with only a few different coherent states. We 
have reduced the problem of Gaussian QPT to Gaus- 
sian QST, for which efficient methods are experimentally 
available [Sg. We have extended our results to multi- 
mode Gaussian QPT and demonstrated that the number 
of test states required increases only polynomially with 
the number of modes. 
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Appendix A: Ancilla-assisted quantum process 
tomography 

Here, we review the existing result of the anciha- 
assisted QPT with TMSV [^ to show some similarities 
but also crucial differences in comparison to our proposal 
of ancilla-free QPT. 

A TMSV is defined by \x{q)) = c, cxp(gat&t)|00), 
where Cq = -^/l — q'^, and q is real. The (unnornialized) 
maximally-entangled state here is 

oo 

|$+) = lim cxp(ga^5t)|00) = V \kk), (Al) 

fe=0 

where a^ (6^) is the creation operator for mode a (b). 
Note that entanglement is not required for the ancilla- 
assisted QPT, but makes it more efficient. In particular, 
the use of the maximally-entangled states can make the 
QPT experimentally optimal with regard to perfect non- 
local correlations Q. 

Assume now that the black box process acts only in 
mode b of the bipartite state \x{q))- After the process, 
we obtain a two- mode state Vlq. One can define the pro- 
jection operator 

T((7) = c,exp[(lnq)a^a], (A2) 

which has the property [40|: 

T{q){a,a^)T-\q) = {a/q,qa^). (A3) 

The TMSV \xiq)) can be written as 

\xiq))=Tiq)®I\^+). (A4) 

According to Eq. (|A4I) . we have 

nq = T{q) ®I-pe- T{q) ® I, (A5) 

where pe = / £ (|$+)($+|). Naturally, 

p,^T-\q)®I-nq-T'\q)®I. (A6) 

We now also formulate the output state of any single- 
mode input state \i^({ck})) = '^k'^k\k) of mode b. Ob- 
viously it can be written as 

= tra(|V*)(V'*l®/-|<J'+)($+l), (A7) 

and \i^*)a = J2k Cfc|fc)a is a single-mode state for mode a 
(sometimes we omit the subscript a or 6 for simplicity). 
We obtain the output state 

Pi, = a(V'*|Pe|V'*)a =tra(|^*)(^*|®I-/3<;) 

= M{{cklq^})\^q\r{{cklq^}))a. (A8) 

More explicit expressions can be obtained by using the 
Husimi Q-function. If the single-mode input state in 



mode & is a coherent state |a), the output state then 
becomes 

p^ ^ (a*|pe|a*) = (a*|T-i(q) ® / • fig • r-i(g) ® I\a*). 

(A9) 
Note that the state \a*) here is a single-mode coherent 
state in mode a. Using the property of T{q) and the 
definition of CSs, a\a*) = a*|a*), we easily find 



T-\q) ® I \a*) ^ Nq{a) \a*/q), 



(AlO) 



where the factor Mq{a) = exp [— |ap(l — l/g^)/2] /cg, 
and \oL* /q) is a coherent state in mode a defined by 
a\a*/q) = {a* / q)\a* / q) . Thus, the output state of mode 
b is 

Pa = |AA,(a)|2 ,{a*/q |fi, \a*/q)a . (All) 

Assume the Husimi Q-function for flq is 
Qn^{Z*,Z^,Za, Zb). According to its definition, 

QnSZlz;,Za,Z^,) = {Za,Zk\nq\Za,Zf,), (A12) 

where \Za,Zi,) is a two-mode coherent state defined by 
(a,b)\Za,Zb) — {Za, Zi,)\Za, Zb). Hence the correspond- 
ing density operator is the following normally-ordered op- 
erator 



n„ 



Qn^ia\bla,b) 



(A13) 



which is simply the operator functional obtained by re- 
placing the variables {Z*,Z^,Za,Zi,) with (0^,6^, a, 6) 
in the Q-function given by Eq. (|A12p . analogously to 
Eq. ^. Therefore, using Eq. (|A11[) and the normally- 
ordered form of fig, wc have the following simple form 
for the Husimi Q-function 

QpjZb*,Zb) = \Nq{a)\^Qn^{a/q,Zi,*,a*/q,Zi,) (A14) 



of the output state pa- Eqs. (jAlll) - (|A14p are the explicit 
expressions of the output state for the input of any co- 
herent state \a). According to Ref. [l^l, if we know the 
output states for all input CSs then we know the output 
states of all states in Fock space. In this approach, given 
any input state \ip), wc can write it in its linear superpo- 
sition form in the coherent-state basis, and then obtain 
the Q-function of its output state by using Eq. (|A14[) . 

These results can be generalized for a multimode 
QPT. To apply the Jamiolkowski isomorphism [41|, we 
consider k pairs of maximally-entangled states, each 
in modes ai,fei, 02, &2,''', o,kTbk- Explicitly, l*!"^) = 
\(j)+)i\(b+)2---\(b+)k. Here |0+), = lim.^i exp((7al6t)|00) 
indicates a maximally-entangled state in modes a,;, bi. 
Subspaces a and b each are now fc-mode. Any state 
■0) in subspace b, can still be written in the form of 
Eq. (jA7p . with the new definitions for \ip) and | $"''}. Us- 
ing Eq. (|A6I) , it is obvious that the output state of these 
fc-pairs of TMSV fully characterizes the process. 
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